Quantum phase transitions and collapse of the Mott gap in the d = 1 + e 
dimensional half-filled Hubbard model 

Jun-ichiro Kishinc* 
Department of Theoretical Studies, Institute for Molecular Science, Okazaki 444-8585, Japan 

(February 1, 2008) 

CT^ , We study the low-energy asymptotics of the half-filled Hubbard model with a 

circular Fermi surface in d = 1 + e continuous dimensions, based on the one-loop 
renormalization-group (RG) method. Peculiarity of the d — 1 + e dimensions is 
incorporated through the mathematical structure of the elementary particle-partcile 
(PP) and particle-hole (PH) loops: infrared logarithmic singularity of the PH loop is 
f~) , smeared for e > 0. The RG flows indicate that a quantum phase transition (QPT) 

(f) ■ from a metallic phase to the Mott insulator phase occurs at a finite on-site Coulomb 

repulsion U for e > 0. We also discuss effects of randomness. 
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I. INTRODUCTION 



Correlation-driven metal-insulator transition (MIT) in the Hubbard model has been a basic prob- 
lem in condensed matter physics JdCentral to this issue is the problem of quantum fluctuation con- 
trolled by the on-site Coulomb repulsion U, the band width W, and the carrier concentration n. In 
the half-filled (n — 1) case, the zero temperature band width control MIT occurs at a critical ratio 
of U = U/W, U c . The exact solution of the Hubbard model is available only in d = 1 dimension, 
where the half-filled gsound state is always an insulator with a finite charge excitation gap (Mott 
gap) at U > U c — 0.0 In d — 2, it is believed that there exists a finite U c except in the case of 
the perfect nesting where U c = 0.B In the case of the d = 2 half- filled Hubbard modpl with the 
nearest neighbor and the second nearest neighbor hoppiags, a Monte Carlo simulation^ a Hartree- 
Fock approximation, and the Gutzwiller approximation^] all indicate that U c is finite. In d = oo, 
the dynamical mean field approachO implies that the quasiparticle spectral weight in the vicinity of 
the Fermi surface vanishes continuously as U approaches a critical value tJ c ~ 1 from belowE The 
0^ | filling control MIT, .-which occurs as the carrier concentration approaches half filling, has also been 

extensively studied .El 

In the d = 1 Hubbard model at half filling, low-energy asymptotics is also well understood in 
terms of the renormalization-group (RG) flow of two-particle scattering strengths, gi,_aa, and 53, 
which correspond to the backward, forward and 2fcF-umklapp scatterings, respectivelyBB The RG 
scheme in d = 1 is based on infrared logarithmic singularities of elementary particle-particle (PP) 
and particle-hole (PH) loops which have the same magnitude and opposite signs. In this context, 
the source of the Mott gap is the umklapp scattering which becomes a relevant perturbation for 
O ■ U > 0. The RG flow also indicates that the charge stiffness reaches zero during the renormalization 

process and the system becomes an insulator. Thus, the RG-based scenario is in perfect agreement 
with the exact solution. 

Dimensionality effects on the d — 1 Mott insulator phase were phenomenologically treated by 
cutting off either of the PP or the PH loop below some characteristic temperatureJij'EJ In the case of 
weakly-coupled chains, dimensional crossovers caused by an interchain hopping t± have been studied 
by treating t± perturbatively and assuming that the scaling procedure in the one-dimaasional regime 
at high energy scales ( uj ^$> t±) remains valid down to the crossover energy scales.E3~til Recently, 
the filling coplrol MIT in coupled Hubbard chains with infinitely large coordination numbers was 
also studied.tj However, these attempts have not clarified the dimensionality effects on the Mott 
gap, because feedback effects of the interchain processes on the Mott gap have been missing. 
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The RG method is straightforwardly extended to the case of d = 1 + e (0 < e <C 1) dimensions .113 
In thiS|-ease, only the PP loop remains logarithmically singular, while the PH loop is smeared for 
e > OO By taking this fact into account from the beginning of renormalization processes, the 
feedback effects on the Mott gap may be incorporated. In this paper, by using the one-loop RG 
method, we examine dimensionality effects on the RG flow of the umklapp process and discuss 
possible QPTs in the Hubbard model with a circular Fermi surface in d = 1 + e dimensions. 

If randomness exists, a QPT between the Mott insulator and a randomness-driven Anderson 
insulator mighj_arise. In the case of d — 1 at half filling, based on the RG method, Fujimoto 
and KawakamO found ithpti sufficiently strong random forward scattering destroys the Mott gap. 
Recently, Ohtsuka et aZ.tSc3 studied the half-filled Hubbard model containing site randomness by 
using the quantum Monte Carlo technique and found that the strong randomness destroys the Mott 
gap. In this paper, we also discuss the QPTs in the half-filled random Hubbard model in d = 1 + e 
dimensions. 

This paper is organized as follows: In Sec. II we introduce the g-ology effective action, derive the 
one-loop RG equations and discuss possible QPTs in the absence of randomness. The effects of 
randomness are discussed in Sec. M, followed by concluding remarks in Sec. IV. 



II. HALF-FILLED HUBBARD MODEL IN D = 1 + e DIMENSIONS 

In this section, we study interplay of electron correlation and dimensionality effects in the half- 
filled Hubbard model with a circular Fermi surface in d = 1 + e dimensions. 



A. Effective action 

We start with the effective action, 

^Hubbard = Y.J ^ J J^G-'iKieKiK^iK) 
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where c* a (K) and c a {K) are the Grassmann variables representing the electron with the spin a 
and K = (k,e) with k and e being a d-dimensional momentum and a Fermion thermal frequency, 
respectively. The non-interacting one-particle propagator is given by 

g- 1 (k,k)=ie-^k). (2) 

The one-particle dispersion is £(fc) = (fc — kp)/2m and the Fermi surface is a d-dimensional sphere 
| k \— kp. Since we consider only the energy scale which is much smaller than the Fermi energy, we 
linearize one-particle dispersion as 

ak)=v F (\k\-k F ), (3) 

where vp is the Fermi velocity. The band width cut off Eq is introduced and the one-particle 
processes are restricted to — Eq < £(k) < Eq. In the two-particle scattering part of ^Hubbard) 
| G\= and | G\= ikp for the normal and the umklapp processes, respectively. It is reasonable to 
assume that the two-particle scattering processes which enter the RG equations in d = 1 + e are, as 
in the case of d = 1, the backward scattering with large momentum transfer | k 3 — k 2 |~ 2kp, the 
forward scattering with small momentum transfer | k$ — k 2 |~ 0, and the 2fc^-umklapp scattering 
characterized by | k^ + k^ — k 2 — k\ |~ 4fc^. We represent the corresponding scattering vertices in 
Figs. l(a)-(c). Dimensionless scattering strengths for these processes are denoted by g%, g 2 , and 33, 



respectively. Unrenormalized scattering strengths are related to the on-site Coulomb repulsion, U , 
as@ 

<?1;0 = .92:0 = .93:0 = U/lTV F = U. (4) 



B. The elementary particle-particle and particle-hole loops 

The peculiarity of d = 1 + e dimensions is incorporated only through the integration measure, 
/ d d k/(2Tr) d . For our purpose here, it is sufficient to integrate over k =| fe| and the angle 9 spanned 
by fc and another fixed momentum. Then we can useEII 

d k , x S d -i f u d-ij U f jQ{ r -„ a \d-2, 



where Sd — 2Tr d / 2 /T(d/2) is the surface area of the d-dimensional unit sphere. 

As is well known, in any dimension, the real part of the elementary PP loop [Fig. 1(d)] at the zero 
total momentum, 

A , , f d d k e(e_ fc )-9(-e fc ) 
J ( 27r ) d w - i_ k - Cfc + 10+ 

with 0(x) being the step function, exhibits an infrared logarithmic singularity of the form, 



d-l 



RAo(w) = - 5d l k ' p log(v/E ). (7) 

(27r) a vp 



In d = 1 + e, we obtain 



sftA (w) ~ --J— ]og(w/Eo), (8) 

Zttvf 

which exactly reproduces the result in d = 1. 

On the other hand, the real part of the elementary PH loop at 2/cf momentum transfer[Fig. 1(e)], 

n , A [ d d k e(^ +Q )-e(c fc ) 
j (2Tr) d uj - i k+Q + e fc + io+ 

with |Q|= 2fc^ , unlike the case of d = 1, no longer exhibits an infrared singularity for <i > 1 and in 
d = 1 + e, takes the form, 
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(10) 



where a) = to /2vpkp and C e is a constant independent of u>. Although this form has already been 
suggested in Ref. [16], we confirm, in the appendix, that the w-dependent term in (Ww is uniquely 
determined. 



C. One- loop renormalization 

One-loop renormalization of the scattering strengths comes from the vertex correction diagrams 
represented in Figs. 2(a) and 2(b) for the normal [g\, g 2 ] and umklapp [g 3 ] processes, respectively. 
The renormalized scattering strength, g[, g' 2 , and g' 3 , are thus given by 



9'i =9i+ .91.92 In — - (.92 - gi)gin(uj), (11) 

.92 = .92 + 2 (.9? + ,9 2 2 ) In |f- - ^(gl + gfMw), (12) 

93 = .93 + 93(91 - 2g 2 )n{uj), (13) 

where 7r(w) = 2irvF$m2k F (u). Fig. 2(a-l) contains the PP loop and all other diagrams contain the 
PH loop. In particular, renormalization of the umklapp process comes from the PH loop only. By 
differentiating equations (0)-(|13) with respect to the scaling parameter 



we obtain the RG equations 



l=]n(E /u), (14) 
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The PH loop gives rise to the smooth cutoff,! 

e/2 



A,= 



d I \ 



= {^£) exp[-d/2] ~ exp[-eJ/2], (18) 



where the ratio of the two cut off energy scales, Eo/2kpVF, is of the order of unity. In the absence 
of the umklapp process, the RG equations obtained here reproduces those in Ref. [16]. 



D. Renormalization-group flow and QPT 

In the d — 1 half-filled Hubbard model, the charge degrees of freedom are governed by the combi- 
nation of (#3, G — gi — 2(72) with the flow lines (G — const.) 2 — g^ 2 = const. For any finite U > 0, they 
are scaled to (g^ — 00, G* = —00), which implies the Mott gap opens due to the relevant umklapp 
scattering. This RG flow also indicates that the charge stiffness K p — ^/(l + G)/(l — G) reaches 
zero during the renormalization process and the system becomes an insulator. To complement this 
scenario, it is useful to map the charge sector of the d = 1 half-filled Hubbard model onto the 
(l+l)-dimensional sine-Gordon model by using the bosonization technique. The sine-Gordon action 
is of the formH 

S SG = Jd 2 r j^( Vr * p (r)) 2 - ^cosh/S^Mr-)]} , (19) 

where r = (x.t) represents (l+l)-dimensional space-time coordinates, ^ P (r) represents a charge 
boson (holon) field, u p is the holon velocity, rand a is a short distance cut off. By applying the RG 
method directory to the sine-Gordon modest 2 ] we obtain the RG equations 

^= 2(1 -K p )g 3 , (20) 



—2 



/ = -^r9 2 ^ 3 , (2D 



where a — L dpp 3 Jo(p) and 2; = Soe - ' is the space-time cutoff. We thus see that 33 becomes 
relevant for the initial condition K p -q < 1 [corresponding to U > for the Hubbard model] and 
then, accordingly, the charge stiffness K p is scaled to zero. 

Now we consider the case of d — 1 + e dimensions. In Figs. 3(a) and (b) are shown the RG flows 
for U = U/ttvf = 0.08 and 0.04, respectively, in d = f.l. In Fig. 3(c), we show the RG flows in 



terms of 173 and G = g\ — 2jj2 for various U in d = 1.1. It is found that there exists a critical value of 
U, U c — 0.0588. For U > U c , the RG flows exhibit runaway trajectories toward (g% = 00, G* = —00) 
[shaded region in Fig. 3(c)], which implies the Mott gap opens at the low-energy limit just as in the 
case of d = 1. The initial values of gz and G at I — correspond to the points on the line, G = —.93, 
represented by a broken line in Fig. 3(c). In the d = 1 half-filled Hubbard model, 53 and G flow 
along this line (denoted by "ID Hubbard") for any U > 0. 

On the other hand, the RG flows approach the fixed points, (g^ = const., G* — 0) for U < U c . 
Marginal behavior of 53 is in accordance with G = g\ — 2#2 — > [see equation (|17|) ] as I — > 00. 
The smooth cutoff, A; = exp[— el/2], in eq. ( p7| ) causes suppression of (73 during the renormalization 
process. This suppression becomes more conspicuous for larger e. The fixed point G* = 0, corre- 
sponding to the non- interacting value of the charge stiffness, K* — 1, implies that the Mott gap 
collapses and the system becomes metallic at the low-energy limit. Thus a QPT from the metallic 
phase to the Mott insulator phase may occur at U = U C . 

Within the RG-based scheme, it remains debatable how the marginal behavior of gs for U > U c 
modifies the ground state property. Regarding this point, recently the density matrix renormalization 
group (DMRC) method was applied to three Hubbard chains coupled via the interchain one-particle 
hopping t±.t3 As a result, it was found that the Mott gap decreases as t± increases. This numerical 
result strongly supprots the RG-based view given here. 

We here give qualitative discussion on the magnitude of the Mott gap. There is no tractable 
method to quantitatively obtain the magnitude of the Mott gap in d > 1. However, the magnitude 
of the Mott gap is qualitatively given by the energy sca^-tjgap = £^e -igap , at which the umklapp 



scattering strength exceeds unity, 53 = 1 [see Fig. 3(a)] .E3Ej To see how w gap reproduces the Mott 



gap, we compare the [/-dependence of the exact Mott gap in d = l,c3 



tt-'U 2 J\ smh[rj/U\ 

with that of w gap in d = 1 . There is arbitrariness in specification of the linearized bandwidth Eq . In 
Fig. 4, we show the case for E Q — OAvf, where cj gap reproduces A cxact (U) well at least for a weak 
U where the weak coupling RG scheme is valid. 

In Fig. 5, we show a low-energy asymptotic phase diagram. We also show uj gajP /Eo as a function 
of d and U. For a fixed dimension, the low-energy asymptotic phase corresponds to a metal and 
the Mott insulator phases for U < U c and U > U c , respectively. The critical value, U c , increases as 
d increases, suggesting a larger U is required for the Mott gap to open as the dimension increases. 
Accordingly, for a fixed U, u gap / E$ decreases with increasing dimensions and disappears at some 
critical dimension, d c [ for example, d c = 1.345 for U = 0.2 ]. 

Apparently Wgap/i^o grows continuously at U = U c for d ~ 1.4. However, the transition is 
always discontinuous by the following reason. As is seen from Fig. 3(c), 33 exceeds unity during the 
renormalization process and gives a finite value of uj gap /E even for U < U c . Nevertheless, uJ gap /E 
has no physical interpretation for U < U c , since the charge stiffness is scaled to the non-interacting 
value. Therefore, uj gap /Eo grows discontinuously at U — U c . The magnitude of discontinuity at the 
transition point becomes more conspicuous as d increases. At present, it is not clear whether the 
discontinuity is an artifact of the one-loop RG method. To elaborate on this point would require 
the two-loop RG analysis, which is too involved a subject to be treated here. 

We here comment on the relevance of the present study to the d — 2 lattice Hubbard model. In 
this paper, we have studied only the case with a circular Fermi surface. Accordingly, our results 
are not smoothly connected to the d = 2 square lattice Hubbard model. In particular, in the 
case of the perfect nesting, the van-Hove singularity at (±7r, ±7r) points gives rise to the "log- 
square" singularity of the elementary PH loop at the momentum transfer Q = (tt,tt). Then the 

Hartree-Fock solution gives the gapQ A ~ ie _27r v*'' 7 , which indicates that the ground state of 
the system is always an insulator for a finite U. This conclusion has been strongly supported by 
numerical studiesoQ However possibility of the MIT at a finite U in the half-filled lattice Hubbard 
model with a Fermi surface of various geometries has been open to question. In the case of the 
d = 2 half-filled Hubbard model with the nearest neighbor afld the second nearest neighbor hopping 
integrals ( t and t' , respectively ), a Monte Carlo simulationja a Hartree-Fock approximation and the 



Gutzwillcr approximation including the antiferromagnetismU give finite critical values, U/t = 2.5, 
2.064, and 3.902, respectively, for t' jt = 0.2. These results are consistent with the present findings 
that, in the case of a circular Fermi surface, the MIT occurs at a finite U for e > 0. 



III. EFFECTS OF RANDOMNESS 

In this section, we study interplay of electron correlation, randomness and dimensionality effects 
in the d — 1 + e dimensional random Hubbard model at half fillingr3 



A. Effective action for quenched randomness 

The action for the scattering processes by the random potentials, 

^random 



X/A/M.)4(vl*M 



(23) 



is added to the Hubbard action (pi). Here, c<r(x,r) = Y2k e X °cr(k,T) with r being an imaginary 
time and v(x) — J^q e iqx v(q) is a random potential at the position x. We assume that the random 

scattering processes which enter the RG equations in d — 1 + e are, as in the case of d = l,c3 
characterized by real and complex random fields r/(x) and £(a:) for the forward scattering with small 
momentum transfer | q |~ and the backward scattering with large momentum transfer | q |~ 2kp, 
respectively, due to the random potential. The random potentials are assumed to be governed by 
Gaussian distributions, 



Prj oc exp 
P$ oc exp 



-D v / d d xn{x) 

d d x£(x)C(x) 



-D 



which lead to 



{r}{x)r]{y)) Ta . nAo 

{i(x)C{y))^ do 



D 



D< 



?S(x-y), 
5(x-y), 



(24) 
(25) 

(26) 

(27) 



where D v = (ttNfTjj)^ 1 and D^ — (nNpT^)^ 1 with r n ^ and Np being the elastic-scattering mean 
free time and the one-particle density of states, respectively. 

We consider the quenched randomness where averaging the free energy is accomplished by means 
of the replica trick which is based on the identity 



?N 



In Z = lim 



N^O N 



(28) 



We introduce N identical replicas of the system labeled by the index a. Then, by using the path- 
integral representation of the partition function, we have 



N 

Z N = j[[ v c a *Vc a exp 



A' 



E*° 



(29) 



where S a = Sji ubbard + "^random ^ s t ne total action and T> symbolizes the measure over the fermionic 
Grassmann variables c a * and c a depending on a replica index a. The replica trick consists of per- 
forming the Gaussian ensemble average (Z N ) lan dom for integer N, continuing the result analytically 
to real N, and taking the limit N — > 0. We thus obtain 



N 



. . - N I" N "I . N N 

drjP v d£dCPt / JJ Vc a *Vc a exp ^S a = J[ Vc a *Vc a exp[^ S a ], (30) 

J J J a —\ \_a=\ J a =\ a = \ 



where the random scattering parts contained in S a are given by 

Jdu r rid 1,1 



n N r 

t^EE/^ 

f3=lo,a' J 


/3=l<r,<r' 



/ dr2 /(^/^ cr(fc ' T1K(fc ' T1)cfr(fe ' r2)cf ' (fc '' T2) 

/ dr2 /(0?/^ cr(fc + O ' ri)c " (fc ' ri)c - r(fe '" O ' T2)c -' (fc '' r2) ' 



where |Q|~ 2fc^. 

We here change imaginary time variables t\ and t-i into Ar = t\ — t 2 and r = {j\ + T2)/2. In the 
integration over At, we introduce a short distance cutoff A and keep only the region, vp | At |< A, 
which couples the two-particle scattering processes and contribute to the RG equations. Then, the 
random forward and backward scattering parts are written as 

5a D V A^^ f d d k f d d k' A [°° d £l 

»n m dom~ -^T^Z^y (2^7 (2^11y 2-f^+S3- £2 -£l) 

c r(fe )£4 )cfr(fe / l£ 3)c^(fe / l£ 2)c«(fe, £ i) (si) 

D.A ^^ [ d d k f d d k' -A- f°° de % .. 

cr(fc + Q>£4)cfr(fc'-Q,e 3 K(fc,e 2 ) C f,(fe', £l ), (32) 



where c"(a:, t) = T 1 / 2 J^ e l( - fe:E er ^c"(fc, e). The actions for the random scatterings inside the same 
replica \J3 — a] are absorbed into the two-particle backward and forward scatterings by introducing^ 

gi=gi-Di, (33) 

52 =52 -A,, (34) 

where D^ — D^A/itvp and D v = D v A/ttv f . Now, in addition to the two-particle scattering vertices 
[Figs. l(a)-(c)], there appear inter-replica vertices as shown in Figs. 6(a) and (b). 



B. One-loop renormalization 

We obtain the vertex correction diagrams for g\, c/2 and g 3 merely by replacing g\ and 172 in Fig. 2 
with g\ and gi-, respectively. However, we must avoid counting the diagram as shown in Fig. 6(c) 
which apparently renormalizes gi, but vanishes in the replica limit, N — > 0, since summation over 
the replica indices 7 = 1,2, ..., N of the inner loop yields N. Keeping this point in mind, we obtain 
the renormalized scattering strength, g[, g' 2 , and g' 3 , which are analogous to equations ([ll|)-(|l3|) and 
given by 

g'i =9i+ 5i52 In ■=- - [(52 - 5i)5i + DI]tt(uj), (35) 

52 = 52 + 2(9? + 92 2 ) In £■ - g (5 2 2 + flf )tt(w), (36) 

53 = 53 + 53(51 - 25 2 )tt(w). (37) 



Renormalization of the inter- replica vertices for f3 ^ a in Figs. 6(a) and (b) comes from the 
vertex correction diagrams as shown in Figs. 7(a) and (b), respectively. We must avoid counting the 
diagram Fig. 7(b-4) which vanishes in the replica limit. We obtain 

D' n = D v -±DJln^-, (38) 

D' s = Ds- AA hi £ + [(251 - &)£* + 2D£]tt(w). (39) 



When we differentiate (p5[)-(p9|) with respect to the scaling parameter, I = \u{Eq/oj), the length-scale, 
A, must also be scaled in accordance with the change of the energy scale, u, as 



dA duo 

— + — = 0. 40 

A uj y ' 



Thus we obtain the RG equations 



*$j. = -£>t- .9132 + [(h - 91)91 + 6|]Ai, (41) 

dff2 
dZ 

^^-g 3 ( gi -2g 2 )X l: (43) 

AD 



D v - (g( + %)/2 + (gi + g 2 3 )Xi/2, (42) 



■n _ n 1 n2 , 



— =D v + Dl/2, (44) 

£) s + ^A, - [(2.91 - 0j)Z) £ + 2D|]A;. (45) 



d£> c 

~dT 



C. Renormalization-group flow and QPT 

i. 7n i/ie absence of the random backward scattering 

First, we consider the case where the random forward scattering is present (-D^ ; o 7^ 0), but the 
random backward scattering is absent (-Df ; o = 0), where D^q and D^-q are initial strengths of the 
random forward and backward scatterings, respectively. In this case, the RG flows indicate that 
the QPT occurs from a metallic fixed point, (D* = 00,53 = 0), to the Mott insulator fixed point, 
(D* = 00, <?3 = 00), as U increases. Typical flows are shown in Figs. 8(a-l) and (a-2) for U — 0.1 and 
0.4, respectively, in the case of d — 1.1 and D v -q = 0.02, where the critical value of U is U c ~ 0.330. 
This behavior in d = 1 + e qualitatively reproduces the case of d = 1,Ej where sufficiently strong 
random forward scattering destroys the Mott gap. 

In Fig. 9(a) is shown a low-energy asymptotic phase diagram, where we also show energy scales of 
the Mott gap, w gap /i?o, introduced in the previous section. We see that both the random forward 
scattering and the raising dimensionality tend to destroy the Mott gap and consequently widen the 
metallic region as compared with the pure case [the phase boundary in the pure case is shown by 
the gray solid line] . 

2. Effects of the random backward scattering 

Next, we consider the case where both the random forward and backward scatterings are present: 
D n -.o 7^ and D^,o ^ 0. The random backward scattering makes it possible for the Anderson 
localization to occur. In this case, there occurs a transition from the Anderson insulator fixed point, 
(Dt = 00, -D* = 00,53 = 0), to the Mott insulator fixed point, (Dt — 00, D* — 00,93 = °°)> as ^ 
increases. Typical flows are shown in Figs. 8(b-l) and (b-3) for U — 0.1 and 0.4, respectively. We 



here introduce the scale l\ oc at which D^ reaches unity [see Fig. 8(b-l)]. Then, oj\ oc — Eoe~ lloc gives 
a qualitative energy scale around which a crossover to the Anderson insulator occurs. In the flows 
of the type of Fig. 8(b-3), 173 always dominates D^ and reaches unity at the scale I = l gap < l\ oc , 
indicating that the Mott gap formation overwhelms the Anderson localization. 

We also find flows toward (£)| = 0, D* = 00. g% = 0), as shown in Fig. 8(b-2). This type of flows 

is found only for < d < 1.575 in the narrow region of U in between the regions corresponding to 
Figs. 8(b-l) and 8(b-3). In these cases, however, D^ exceeds unity at some scaling parameter Zi oc 
during the renormalization, which indicates that the perturbative treatment breaks down and the 
localization occurs around the energy scale specified by I = l\ oc . Thus we interpret that the ground 
state corresponding to this fixed point is the Anderson insulator. It is beyond the RG-based scheme 
to settle this ambiguity and we do not go into the details on this issue here. 

In Fig. 9(b) is shown a low-energy asymptotic phase diagram, where we also show wi oc / Eq and 
<x> g ap/-Eo- As compared with Fig. 9(a), the phase boundary remains nearly unchanged, but the 
metallic phase in Fig. 9(a) is replaced with the Anderson insulator phase due to the random backward 
scattering. 

The present results indicate that the QPT from the Apderson to the Mott insulators occurs in 
both d = 1 and d > 1. Recently, Ohtsuka and Hatsugaic9 studied the half- filled Hubbard model 
containing site randomness by using the Monte Carlo method. They found that the QPTs from an 
incompressible (Mott) to a compressible (Anderson) insulators occur in all the cases of d — 1,2,3. 
This numerical result is consistent with the RG-based views given here. 



IV. CONCLUDING REMARKS 

In this paper, based on the one-loop RG flows, we have studied QPTs in the half-filled Hubbard 
model with a circular Fermi surface in d — 1 + e continuous dimensions. Peculiarity of the d = 1 + e 
dimensions was incorporated only through the mathematical structure of the elementary PP and 
PH loops: infrared logarithmic singularity of the PH loop is smeared. We have studied the following 
three cases: 

(1) In the absence of randomness: The QPT from the metallic phase to the Mott insulator phase 
occurs at a finite U for e > 0. 

(2) In the case where the random forward scattering is present, but the random backward scattering 
is absent: Both random forward scattering and raising dimensionality tend to destroy the Mott gap. 
Consequently, U c becomes finite for e > and the metallic region becomes wider as compared with 
the pure case (1). 

(3) In the case where both the random forward and backward scatterings are present: The phase 
boundary remains nearly unchanged as compared with the case (2), but the metallic phase in the 
case (2) is replaced with the Anderson insulator phase due to the random backward scattering. 

In the present study, the ground state properties were conjectured based solely on the orierloop 
RG flows. At present, numerical studies are in progress to complement the views given hereB3 
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APPENDIX A: DERIVATION OF (10) 

Our purpose here is to show that the w-dependence of 5Rn 2 fe F (w) is uniquely determined as (|l0[) . 
We start with equation (pi). The imaginary counterpart is given by 



^U 2k Au J )^~n^Jk d - 1 dkj 7r d9( S in6) d - 2 e(-^ k )e^ k+Q )S(uj-^ k+Q +^ k ), (Al) 

for < u) < 2vpkp, 3n2fc F (w) = for 2vpkp < w and satisfies 3ir 2 fc F (— u>) = — 9n 2 fc F (cj). The 
delta function in the integrand of (Al) is rewritten as S(u> — £fc+a ^ £fc) ~ v\ ( \ + T ) ^ ~ *°)> 



where we introduced Q = to/2vFkp, k = k/kp, t = cos9 and to = w — (1 — Co )/k. So we have 



(2tt)<V J^z \2 k 

= -*t§C a - ^ r dn ( " 2 - **> <" 2 - ^ <A2> 

where u = fc + u). For small w, since the region of integration is limited to the vicinity of u = 1, it 
is reasonable to replace the integrand, (u 2 — Co 2 ) (u 2 — I) -5- , with 2 _ 2~(w — l) -5- and we have 

-""IsT^^ 1 -* 2 ^*^ < A3 > 

2(27r) a Ui? a— 1 

which corresponds to a special case of equation (7.12) in Ref. [17] for | q |= 2kp. In rf = 1 + e, noting 
Sd-i ~ e, we obtain 

3h 2fcf ( w ) ~ - J-(l - &)-W&/*. (A4) 

4wf 

The real counterpart is obtained by Kramers-Kronig transformation, 

smi^H^^/ 1 ^ "?; 2 ^ , (as) 

where the symbol V denotes Cauchy pri ncip al value integral [Note that < w < 1]. Since the 
prefactor, (1 — a) 2 )~ 1+e / 2 , on the r.h.s o f (lA4|) has already appeared in the exact expression, (|A2|), 
and the remaining part of the integral, (|A2|), shows no singularity at u> = 1, the expression for small 



Co, ( A4), holds analytical property of 3Il2fc F (to) correctly even for tD~l. Thus it is reasonable to use 



the expression (A4) in (A.5) and we have 



m ^) ~ -^j^ \^.X(i-W^ (A6) 



which is to be evaluated for small Co and e. 



To evaluate (A6), let 



/(*) = ,2 -2U2 tJ ^ - ^^ ( A? ) 

(z 2 — w^)(z 2 — 1) 

and consider the integral, § c f(z)dz, along the contour as depicted in Fig. 10. f(z) has poles 
at z = ±u) and branch points at z = 0, 00 and z = ±1. We choose branch cuts in the region 
{3iz < —1} U {0 < Sftz}. In the limit as the large circle recedes to infinity, it gives no contribution. 
The residue at the pole z — —oj gives 

f(z)dz = -nie™' 2 "f /2 ~ -™ e / 2 . (A8) 

c (1 -w 2 ) 1 £ / 2 

The remainder of the contour is deformed into an integral enclosing the cuts and encircling the pole 
z = Co and the branch points at z = and z = ±1. The points z = 0, ±1 give no contribution. The 
integrals encircling the pole z = Co and the remainder along the real axis give 

10 



f(z)dz = %^ + e^n "*'* + 27rv F (e^ - e^)m 2kF (^) 

c z I 1 ~ ^ ) ' 

, 4„ ie . f°° x[x(x 2 - 1)] £ / 2 , . 3.. »., /- 1 x[~a;(.T 2 - 1)] £ / 2 , 

y x (x 2 - w 2 )(a; 2 - 1) u /_ 00 (a; 2 - w 2 )(.x 2 - 1) 

Here the integrals in the second line are evaluated for small Gj and e as 

x[x{x 2 ~ l)] e / 2 /- 1 x[-a:(x 2 - 1)] £ / 2 



1 (X 2 ~ CJ 2 ){X 2 - 1) /_„,, (.T 2 - £ 2 )(.T 2 - 1) 

, L l n 2r(l-e/2) ' ' 

where the integral in the second line converges for < e < 4/3. We thus obtain, for small Cj and e, 

i /0)dz = 7rid) £/2 - 27r 2 w F ie5Rn 2fcF (w) + 2?ri - 7ri. (A9) 

Therefore we obtain "$X\ik F (oj) of the form which is correct up to the leading order of u>, 

r^l 2 



3f?n 2fcF ( w ) 



2ttvf 



6/2 



a 



(A10) 



which is (p. pp. Although above manipulation gives C e = 1/e, it seems feasible to choose C c = — 2/e as 
suggested in Ref. [16] to reproduce correctly the limit form of $Ol2k F (w) = » ^ logo) at e = 0. This 



discrepancy may arise, because to evaluate (|A5[) we used the expression (|AJ) which holds analytical 
property of QTl2fc F (w) correctly but misses contribution from u)~l. We do not go into the details 
here, since an explicit form of C e does not enter the RG equations (|l5|)-(|l7|) and (|4l|)-(|4q). 
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FIG. 1. Two-particle scattering vertices for the (a) backward, (b) forward, and (c) 2fc_F-umklapp scat- 
tering processes, and the elementary (d) particle-particle (PP) and (e) particle-hole (PH) loops. 



FIG. 2. Vertex correction diagrams for the (a) normal [gi, 52] and (b) umklapp [gs] processes. 
White and black circles represent the normal and umklapp scatterings, respectively. 



FIG. 3. RG flows in the case of d = 1.1 with (a) U = 0.08 and (b) 0.04. (c) The RG trajectories in terms 
of <?3 and G — gi — 2g2 in d — 1.1. A critical value is U c — 0.0588. In the d = 1 half-filled Hubbard model, 
<73 and G flow along the broken line denoted by "ID Hubbard". 



FIG. 4. [/-dependence of A exa ct(t/) and w gap = Eoe igap in d — 1. 



FIG. 5. A low-energy asymptotic phase diagram of the d = 1 + e dimensional Hubbard model at half 
filling. We also show energy scales of the Mott gap, ui ga , p /Eo, as a function of d and U. 



FIG. 6. Inter-replica vertices originating from the random (a) forward and (b) backward scatterings, (c) 
The diagram which is proportional to the number of replicas, N, and vanishes in the replica limit, N — ► 0. 



FIG. 7. The vertex correction diagrams for the inter-replica (a) forward [D v ] and (b) backward [D%] 
processes. We must avoid counting the diagram (b-4) which vanishes in the replica limit. 



FIG. 8. (a) RG flows of D v and g 3 for (a-1) U = 0.1 and (a-2) U = 0.4 in the case of d = 1.1, D vi0 = 0.02 
and Dt-.Q = 0. (b) RG flows of D v , D i and g :i for (b-1) U = 0.1, (b-2) U = 0.2 and (b-3) U = 0.4 in the case 
of d = 1.1, A, ; o = 0.02 and Z) e;0 = 0.08. 



FIG. 9. Low-energy asymptotic phase diagrams of the d = 1 + e dimensional random Hubbard model 
at half filling in the cases where (a) the random forward scattering is present, but the random backward 
scattering is absent, and (b) both the random forward and backward scatterings are present. In (a), the 
phase boundary in the pure case [see Fig. 5] is shown by the gray solid line. We also show uj\ oc /Eq and 
w g ap/i?o, as a function of d and U. 



FIG. 10. Contour to evaluate the integral (A6). 
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